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The  paper  presents  a  new  approach  for  state  estimation  of  lithium-iron  phosphate  batteries.  Lithium— iron 
phosphate/graphite  batteries  are  very  intricate  in  state  of  charge  estimation  since  the  open  circuit  voltage 
characteristic  is  flat  and  ambiguous.  The  characteristic  is  ambiguous  because  open  circuit  voltages  are 
different  if  one  charges  or  discharges  the  battery.  These  properties  also  hinder  state  of  health  estimation. 
Therefore  conventional  approaches  like  Kalman  filtering  which  represents  a  state  by  only  the  mean  and  the 
variance  of  a  Gaussian  probability  density  function  tend  to  fail.  The  particle  filter  presented  here  overcomes 
the  problem  by  using  Monte  Carlo  sampling  methods  which  are  able  to  represent  any  probability  density 
function.  The  ambiguities  can  be  modelled  stochastically  and  complex  models  dealing  with  hysteresis  can 
be  avoided.  The  state  of  health  estimation  employs  the  same  framework  and  takes  the  estimated  state  of 
charge  as  input  for  estimating  the  battery’s  state  of  health.  The  filter  was  developed  for  A123  lithium-iron 
phosphate  batteries.  For  validation  purposes  user  profiles  for  batteries  in  different  ageing  states  like 
electric  vehicles  and  off-grid  power  supply  applications  were  generated  at  a  battery  testing  system.  The 
results  show  good  accuracy. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium-ion  batteries  are  the  most  favoured  battery  technology 
in  many  upcoming  applications  today,  may  it  be  for  electric  vehicles 
or  for  storing  renewable  energy.  Classical  lithium-ion  technology 
though  suffers  from  safety  problems  and  cost  of  raw  materials  as 
e.g.  cobalt  [1],  Therefore  researchers  are  working  on  alternative 
materials  like  phosphate  based  batteries  [2],  This  is  a  very  inter¬ 
esting  class  of  cathode  materials  especially  for  safety  reasons.  But 
phosphate  materials  tend  to  have  very  flat  open  circuit  voltages 
with  a  hysteresis  between  charging  and  discharging.  For  LiFePCU, 
the  most  widely  spread  of  these  materials,  this  characteristic  is 
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shown  in  Fig.  1.  During  partial  cycling  even  micro  hysteresis  appear 
like  one  can  also  see  for  nickel  based  batteries  [3,4].  This  is 
a  problem  for  system  engineering  which  needs  to  estimate  the  state 
of  charge  based  on  the  terminal  voltage.  This  paper  will  introduce 
a  framework  called  particle  filter  which  overcomes  this  problem  by 
modelling  this  behaviour  stochastically. 

First  the  method  will  be  introduced,  afterwards  the  method  will 
be  proven  based  on  real  measurement  data  from  the  batteries.  In 
the  conclusions  the  results  are  assessed  and  further  works  will  be 
described.  For  this  paper  2.3  Ah  batteries  from  A123  are  used 
having  the  type  designation  ANR26650M1. 

2.  Description  of  the  parallel  particle  filter 

The  most  popular  filter  within  the  family  of  Bayesian  filters  is 
the  Kalman  filter  [5—12].  The  Kalman  filter  is  an  analytical  solution 


706 


S.  Schwunk  et  al.  /  Journal  of  Power  Sources  239  (2013)  705-710 


State  of  charge 

Discharge  Charge  - 

Fig.  1.  Open  circuit  voltage  curves  of  a  lithium-iron  phosphate  battery  during 
charging  and  discharging  of  the  batteiy.  One  can  clearly  see  the  very  flat  characteristic 
and  the  hysteresis  between  the  two  curves. 

of  the  Bayesian  filter  for  Gaussian  distributions.  By  employing 
Monte  Carlo  sampling  techniques  the  particle  filter  offers  the 
possibility  to  deal  with  any  kind  of  distribution  by  approximating 
the  respective  probability  density  function  by  a  set  of  particles  or 
samples.  For  LiFePC^-batteries  with  their  hysteresis  in  the  open 
circuit  voltage  this  offers  the  possibility  to  use  multimodal  distri¬ 
butions  modelling  this  behaviour  stochastically. 

2.1.  Recursive  Bayesian  filtering 

A  recursive  Bayesian  filter  estimates  the  state  of  dynamic 
systems  online.  It  takes  into  account  the  probability  for  reaching 
a  state  xt  with  respect  to  the  states  estimated  in  time  steps  xi;t_i 
before  and  the  probability  for  observing  a  certain  quantity  zt  when 
being  in  a  certain  state  xt.  The  system  is  described  by  three 
quantities: 

1.  State  xt :  This  is  the  state  of  the  system  at  time  t,  which  shall  be 
estimated.  The  state  can  neither  be  observed  nor  measured. 

2.  Input  ut:  The  input  will  change  the  system’s  state  over  time  and 
can  be  measured. 

3.  Output  zt:  The  output  is  in  some  correlation  with  the  system’s 
state  xt  enabling  at  least  a  rough  estimation  of  the  state.  The 
quantity  can  be  observed  and  measured. 

The  Bayesian  filter  estimates  the  following  probability  density 
function  based  on  the  quantities  explained  before: 

P(xt|Zi:t,Ui;t_i)  =  (1) 

As  one  applies  the  Bayes’  theorem  [13]  one  gets  the  following 
equation.  The  denominator  normalising  the  function  to  one  is 
written  as  p. 

(  °  P(zt|z1:t  l7u1:t  0  (  } 

=  r/-1P(zt|xt,z1:t-1,Ui:t_i)P(Xt|zi:t_i,Ui:t_1)  (3) 

Now  the  equation  is  marginalised  over  xt_i  and  one  assumes 
a  Markov  chain  [13].  Therefore  only  the  quantities  ut_i  and  zt  and 
not  the  former  values  have  an  effect  on  the  calculated  probability  of 
the  current  state: 


P(xt)  =  V  i P(Zt|Xt)  ■  J  P(xt|X(_i  ,Ut—\)  ■P(Xf_1)dXf_1  (4) 

The  probability  density  function  P(xt_ i)  describes  the  state’s 
probability  density  function  in  the  time  step  before.  The  function 
P(xt|xt_1,  u£_i)  takes  the  influence  of  the  inputs  ut_i  for  the 
progression  of  the  system  from  the  state  xt_£  to  the  state  xt  into 
account.  The  function  P(zt|xt)  represents  the  probability  for 
observing  the  measurement  zt  given  a  certain  state  xt. 

2.2.  The  particle  filter 

Since  for  most  probability  density  functions  there  is  no  easy 
solution  of  this  equation,  techniques  were  developed  for  coping 
with  more  complex  distributions  than  the  Gaussian  distributions  of 
the  Kalman  filter.  In  Refs.  [13,14]  methods  are  introduced,  which 
employ  Monte  Carlo  methods  for  representing  the  distributions. 
The  state  distribution  is  represented  by  a  set  of  samples,  therefore 
being  able  to  represent  any  distribution  though  of  course 
with  accuracy  limited  by  the  number  of  particles  employed.  High 
probabilities  are  represented  by  a  huge  number  of  particles  in 
a  certain  area,  low  probabilities  by  a  low  number  or  no  particles  in 
an  area. 

The  algorithm  is  performed  in  three  steps.  The  first  step  is  the 
state  transition.  In  this  step  the  influence  of  the  input  ut- 1  on  each 
sample  is  calculated.  The  uncertainty  of  this  step  and  the 
measurement  errors  of  the  input  are  taken  into  account  by  the 
addition  of  some  noise  on  each  sample  during  state  transition. 
The  result  is  the  probability  density  function  P(x t|xt_i,  ut_i).  Over 
time  there  will  be  a  diffusion  of  the  samples  increasing  the  variance 
of  the  estimation. 

The  second  and  the  third  step  assess  the  plausibility  of  a  state  xt 
when  observing  a  measurement  zt.  This  leads  to  decreasing  the 
variance  and  limiting  the  diffusion  of  the  particles. 

In  the  second  step  the  samples  are  weighted.  According  to  the 
observed  measurement  value  zt  a  weight  of  each  sample  sf  is 
calculated  via  a  probability  density  function  P(zt |xt).  Afterwards  the 
overall  weight  of  all  samples  is  normalised  to  one. 

The  third  step  is  resampling  all  the  samples  according  to  their 
weights.  After  this  step  all  samples  have  the  same  weight  again.  In 
Ref.  [13]  several  sampling  methods  are  introduced  that  might  be 
employed.  For  low  calculation  effort  low  variance  resampling  is 
employed  [15], 

2.3.  Application  of  the  particle  filter 

In  this  paper  the  particle  filter  is  used  for  estimating  the  two 
most  important  quantities:  the  state  of  charge  of  the  battery  and 
the  state  of  health  of  the  battery.  The  state  of  charge  is  defined  as 
the  amount  of  charge  in  the  battery  divided  by  its  capacity: 

SOC  =  J1  (5) 

'-batt 

The  state  of  health  of  the  battery  is  a  normalised  value  defined 
as  the  battery’s  actual  capacity  Cbatt  divided  by  its  nominal 
capacity  Cr: 

SOH  =  (6) 

Gr 

Each  of  the  quantities  is  represented  by  a  set  of  particles  of  size 
^samples*  The  sizes  of  the  sets  may  differ.  For  both  quantities 
a  separate  particle  filter  is  implemented  and  they  exchange  the 
averages  of  the  samples  representing  the  quantity.  The  state  of 
charge  filter  takes  the  estimated  state  of  health  as  a  parameter  to  its 
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model  during  state  transition,  which  is  only  updated  in  long  time 
intervals.  The  state  of  health  filter  takes  the  state  of  charge  as  the 
output  zt  of  the  system,  which  can  be  observed  for  weighting  the 
samples,  which  represent  the  state  of  health. 
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^discharge, t  ~  O^discharge  (sS0C,t)  +  (sSOC,t  >  4>att)  '  4>att  ( 1 1 ) 
^charged  =  0CVcharge  (sSOC,t)  +  Ri  (sSOC,t>  4>att)  '  4>att  ( 1 2) 


2.4.  Estimating  the  state  of  charge 

The  state  of  charge  is  the  state  xt  of  the  system.  Before  the 
particle  filter  starts  working  all  samples  need  to  be  initialized. 
Because  the  actual  state  of  charge  is  generally  unknown,  all 
samples  are  equally  distributed  over  the  complete  state  of  charge 
range.  During  state  transition  the  following  equation  calculates 
the  influence  of  the  battery  current  /batt  under  consideration  of 
a  noise  z  representing  the  modelling  and  measurement 
uncertainties: 


Both  voltages  are  compared  to  the  measured  terminal  voltage 
and  two  probabilities  are  determined  by  the  Cauchy  distribution 
and  added  to  form  the  weight  of  each  sample: 
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The  noise  r  is  sampled  from  a  suitable  distribution.  For  avoiding 
too  little  diffusion  a  Cauchy  Lorentz  distribution,  which  has  very  fat 
tails  compared  to  a  Gaussian  distribution,  is  employed: 
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1  s 

tt  s2  +  (x  —  x0)2 


(8) 


For  sampling  from  this  distribution  one  draws  samples  c 
uniformly  [0,1  ]  and  solves  the  following  integral  for  receiving  the 
nonuniformly  distributed  samples  y  [13]: 


c  = 


(9) 


For  a  Cauchy  Lorentz  distribution  the  following  equation 
transforms  the  uniformly  distributed  into  Cauchy  distributed 
samples: 

y  =  x0  +  s-tan(7r- (c  -  0.5))  (10) 


The  addition  of  noise  leads  to  a  diffusion  of  the  samples  making 
the  estimation  more  inaccurate  over  time.  For  decreasing  the 
diffusion  again  a  model  is  employed  making  it  possible  to  compare 
the  samples  to  a  measured  value  and  checking  each  of  them  for 
plausibility.  In  case  of  the  state  of  charge  estimation  a  quasista¬ 
tionary  battery  model  is  employed.  Since  this  is  a  LiFeP04  batter- 
y  the  hysteresis  makes  the  possible  state  of  charge  values 
ambiguous,  resulting  in  several  possibilities  for  the  correct  state  of 
charge  value. 

This  is  stochastically  modelled  by  the  function  P[zt\xt),  in  this 
case  a  bimodal  Cauchy  distribution.  Two  possibilities  for  the  ex¬ 
pected  terminal  voltage  for  each  sample,  I'discharge  and  Vcharge  are 
calculated  by  the  following  equations  for  the  discharge  and  the 
charge  case  of  the  open  circuit  voltage: 


The  long  tails  avoid  sample  impoverishment  and  lead  to  higher 
stability  and  robustness  of  the  algorithm.  This  results  in  a  set  of 
couples  consisting  of  a  sample  and  its  probability.  These  probabil¬ 
ities  are  normalised  to  one  and  then  resampled  using  the  low 
variance  resampling. 

The  low  variance  resampling  puts  all  the  samples  in  one  line 
assigning  each  sample  a  length  according  to  its  weight.  For  starting 
the  resampling  process  a  random  starting  value  is  chosen  within 
the  interval  [0,1  /N)  in  which  N  is  the  number  of  samples.  Afterwards 
the  algorithm  goes  forward  with  step  size  1  /N  in  N  -  1  steps  and 
chooses  at  each  step  the  corresponding  sample.  Thereby  exactly  the 
same  amount  of  new  samples  as  the  old  sample  set  is  generated 
and  samples  with  small  weight  are  neglected  and  samples  with 
a  high  weight  are  chosen  several  times.  The  algorithm  is  illustrated 
in  Fig.  2. 

Since  a  bimodal  Cauchy  distribution  is  used  also  the  Monte  Carlo 
representation  can  be  multimodal  as  is  illustrated  in  Fig.  3.  For  the 
state  estimation  the  mean  of  the  particles  is  calculated  and  used  as 
estimation. 

2.5.  Implementation  of  the  state  of  health  filter 

The  state  of  health  is  the  state  of  the  system  xt.  The  state  of 
health  filter  assumes  no  inputs  that  influence  the  capacity  of  the 
battery,  resulting  in  P(xt|x£_i,Ut_i)  =  P(x£|x£_!).  A  very  slow  diffusion 
of  the  particles  is  assumed  for  modelling  the  capacity  change  over 
time.  For  the  diffusion  a  noise  drawn  from  a  Cauchy  distribution  is 
added  to  each  particle: 

sSOH,t  =  sSOH,t-l  +  ^  (14) 

For  reacting  on  the  diffusion  processes  the  output  of  the 
state  of  charge  filter  is  used  to  calculate  the  measurement  value 
ASOCmeasured  and  compares  it  to  the  integrated  charge  flown  into 
and  out  of  the  battery  Qstep-  Therefore  the  difference  of  state  of 
charge  between  two  points  in  time  is  the  measurement  value: 
zt  =  ASOCt.  Each  of  the  particle  is  weighted  by  a  Cauchy  distribution 
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Fig.  2.  Low  variance  resampling  for  generating  a  new  unweighted  sample  set  from  a  weighted  sample  set. 
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State  Of  Charge:  0.487441  /Time:  65.5259  h 


State  of  Charge 


Fig.  3.  Histogram  of  the  particles  during  cycling  in  a  medium  and  ambiguous  state  of 
charge.  One  can  clearly  see,  how  a  bimodal  distribution  evolves  due  to  the  modelled 
hysteresis  of  the  open  circuit  voltage. 


and  assigned  a  weight  accordingly.  The  expected  state  of  charge 
change  is  calculated  for  each  particle  describing  the  state: 

ASOCf  =  ^  (15) 

SSOHt'<-r 


The  weight  for  each  of  these  particles  is  then  assigned  by 
a  Cauchy  distribution: 
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After  assigning  a  weight  to  each  particle  the  resampling  step  is 
performed  with  the  low  variance  resampling  method  and  a  new 
unweighted  sample  set  is  generated.  For  estimating  the  state  of 
health  the  mean  of  the  particles  is  generated. 


Current  - 

Fig.  4.  Current  profiles  for  a  photovoltaic  off-grid  power  supply  and  an  electric  vehicle 
application  as  applied  here  for  validating  the  state  of  charge  and  state  of  health 
estimations. 


3.  Validation 


3.1.  Measurements 

For  validating  the  algorithm  current  profiles  were  run  on 
a  battery  cycler  which  shall  represent  specific  applications.  Two 
applications  were  chosen.  First  a  profile  for  a  photovoltaic  off  grid 
power  supply  is  generated  with  low  currents  but  only  little  breaks 
in  the  consumption.  Second  a  profile  is  generated  with  more  breaks  o 
but  higher  charge  and  discharge  currents,  which  are  based  on 
standard  driving  cycles,  representing  an  electric  vehicle  applica¬ 
tion.  The  photovoltaic  profile  was  performed  on  a  new  battery 
having  the  state  of  health  1.0  and  the  electric  vehicle  profile  on  an 
already  aged  battery  having  the  state  of  health  0.939.  The  ageing 
states  were  determined  in  capacity  tests  before  profiles  were  run. 

The  current  profiles  can  be  seen  in  Fig.  4. 

3.2.  Photovoltaic  profile 

In  Fig.  5  you  see  the  reference  state  of  charge  and  state  of  health 
for  a  lithium— iron  phosphate  battery  with  state  of  health  100%.  At 
the  beginning  the  sample  set  for  the  state  of  charge  estimation  is 
equally  distributed  in  the  interval  [0,1  ]  resulting  in  starting  value  of 
about  0.5  for  the  state  of  charge  estimation.  The  sample  set  for  the 
state  of  health  estimation  is  equally  distributed  in  the  interval 
[0.35,1.35]  resulting  in  a  starting  value  of  about  0.85  for  the  state  of 


SOC  estimation  particle  filter  - 

SOC  reference 

SOH  estimation  particle  filter  . 

SOH  reference 

Fig.  5.  Results  for  a  current  profile  in  a  photovoltaic  off-grid  power  supply.  The  filter 
finds  both  state  of  charge  and  state  of  health  correctly.  The  battery  has  got  a  state  of 
health  of  100%.  The  reference  for  the  state  of  charge  was  generated  by  fully  charging 
the  battery  for  knowing  the  exact  state  of  charge  at  the  beginning  of  the  test  and 
measuring  and  integrating  the  current  afterwards.  For  the  state  of  health  a  capacity 
test  was  performed  before  applying  the  off  grid  power  supply  profile  on  the  battery. 
The  references  are  therefore  both  measured  data. 
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health  estimation.  The  same  value  is  used  as  starting  value  for  the 
state  of  health  used  in  the  state  of  charge  estimation. 

As  can  be  seen  the  state  of  health  estimation  only  starts  after  the 
state  of  charge  has  found  a  reasonable  value.  For  achieving  that 
a  starting  time  for  the  state  of  health  filter  is  introduced.  The  state 
of  charge  filter  quickly  finds  the  correct  value  and  is  able  to  track  it  o 

over  the  whole  time,  while  still  doing  some  corrections  as  can  be  8 

seen  during  the  rest  phases.  This  guarantees,  that  the  state  of 
charge  estimation  is  no  pure  Ah  counter. 

The  filter  tracks  both  estimated  states  very  well  over  time,  is 
stable  and  there  is  hardly  no  noise  on  the  quantities.  This  proves 
that  the  stochastic  modelling  of  the  hysteresis  is  sufficient  to  cope 
with  the  ambiguities  of  the  hysteresis  and  the  resulting  uncer¬ 
tainties  for  the  state  of  charge  estimation. 

For  proving  that  the  algorithm  is  not  in  need  of  starting  with 
a  full  charge,  in  another  run  the  algorithm  is  started  later  at  a  lower 
state  of  charge.  As  can  be  seen  clearly  in  Fig.  6,  the  state  of  charge  is 
found  quickly  and  tracked  well  over  time. 

3.3.  Electric  vehicle  profile 

The  filter  is  parameterised  exactly  the  same  way  for  the  electric 
vehicle  profile  as  for  the  photovoltaic  profile.  An  aged  battery  with 
state  of  health  93.9%  is  used  and  for  this  aged  battery  results  are 
also  very  good  and  robust,  though  the  underlying  battery  model  is 
only  quasistationary  and  developed  for  a  new  battery.  Therefore  the 
ageing  of  the  battery  resulting  in  an  increase  in  inner  resistance  is 
not  captured  by  the  model. 

The  results  for  state  of  charge  and  state  of  health  from  the 
validation  are  shown  in  Fig.  7.  The  state  of  charge  is  tracked  very 
well  over  the  whole  time  providing  a  good  basis  for  the  state  of 
health  estimation.  Therefore  also  the  correct  state  of  health  is  found 
quickly  and  then  stably  tracked  over  time. 

3.4.  Non-deterministic  behaviour 

In  Fig.  8  the  non-deterministic  behaviour  of  the  particle  filter  is 
shown.  As  can  be  seen  each  state  of  charge  estimation  and  state  of 
health  estimation  is  different,  though  the  data,  the  algorithm  and 
the  parameters  are  the  same.  The  only  difference  is  the  seed  for  the 
function  generating  the  equally  distributed  random  number.  Those 
are  used  for  the  noise  added  during  the  calculation  the  effect  of  the 


SOC  estimation  particle  filter  - 

SOC  reference 

SOH  estimation  particle  filter  . 

SOH  reference  . 

Fig.  7.  Results  for  applying  the  particle  filter  to  an  electric  vehicle  profile.  The  battery 
has  only  got  a  state  of  health  of  93.9%.  The  reference  for  the  state  of  charge  was 
generated  by  completely  charging  the  battery  to  know  the  exact  state  of  charge  at  the 
beginning  of  the  test  and  measuring  and  integrating  the  current  afterwards.  For  the 
state  of  health  a  capacity  test  was  performed  before  applying  the  electric  vehicle 
profile  on  the  battery.  Both  references  are  therefore  measured  data. 
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Fig.  6.  The  particle  filter  is  started  later  on  the  same  data  set  as  for  Fig.  5,  starting 
therefore  only  at  a  state  of  charge  of  80%.  As  can  be  seen  clearly  no  full  charge  is 
needed  for  a  good  state  estimation. 


Fig.  8.  Different  runs  of  the  identical  particle  filter  for  the  same  data,  in  this  case  for  an 
off-grid  power  supply.  As  can  be  seen  clearly  the  result  differs  from  run  to  run  due  to 
the  random  numbers  drawn  by  the  Monte  Carlo  method. 
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inputs  on  the  states  and  for  choosing  the  starting  value  for  the  low 
variance  resampling. 

4.  Conclusions 

A  framework  for  dealing  with  difficult  and  ambiguous  batteries 
like  LiFePCU/graphite  batteries  was  presented.  The  ambiguous 
range  of  the  open  circuit  voltage  is  stochastically  modelled  opening 
up  possibilities  of  multimodal  estimations  for  the  battery’s  state  of 
charge  by  a  Monte  Carlo  representation  of  its  state  distribution  with 
samples.  Thereby  the  algorithm  is  both  able  to  find  a  correct  value 
at  the  start  and  track  the  state  of  charge  correctly.  This  forms  the 
basis  for  the  good  state  of  health  estimation. 

The  main  advantage  of  this  stochastic  modelling  compared  to 
a  physical  approach  for  tracking  the  hysteresis  behaviour  is  a  quick 
modelling  during  development  time,  no  starting  state  for  the 
hysteresis  needed  and  a  possibility  for  indicating  several  plausible 
states  with  their  probabilities.  The  last  issue  is  especially  inter¬ 
esting  during  the  start  up  phases. 

The  algorithm  reaches  a  high  accuracy  in  the  examples  shown  in 
the  paper.  The  error  values  for  state  of  charge  and  state  of  health 
estimation  settle  at  very  low  levels. 

Future  work  will  focus  on  analysing  and  increasing  robustness  of 
the  method,  on  extracting  more  information  out  of  the  sample  sets 
like  error  boundaries  and  estimation  of  correctness  and  analysing 
the  long  term  behaviour  of  the  algorithm  during  ageing. 


For  increasing  accuracy  and  accelerating  the  determination  of 
the  correct  state  of  charge  one  could  use  more  complex  impedance 
models  instead  of  the  more  simplified  battery  model  used  in  this 
paper.  In  general  this  algorithm  could  be  a  possible  replacement  for 
the  different  kinds  of  Kalman  filters  used  today  as  state  of  the  art  for 
state  estimation. 
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